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tion of the natural channel radius is introduced. As a result, it is shown to 
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be possible to use well-developed technique of R- and P-matrix theory for 
calculation of resonant states characteristics, scattering phase shifts, etc., in 
the approaches based on harmonic oscillator expansions, e.g., in nuclear shell- 
model calculations. P-matrix is used also for formulation of the method of 
treating Coulomb asymptotics in the scattering theory in oscillator represen- 
tation. 

I. INTRODUCTION. 

A number of methods for treating effects of the continuum spectrum in nuclear structure 
calculations have been developed. We shall mention the i?-matrix approach ||l| (see also 
review paper and references therein) and related to it the P-matrix approach (see IQ and 
references therein), the Feshbach projection method 0] and, based on that, the continuum 
shell model 0,3, early coupled channel calculations |'^^T0| and more developed various 



continuum RPA approaches |Tl|-|13|, methods based on the expansions of the continuum 



spectrum wave function in Sturmian (Weinberg) functions [p!^,p!5| or in other basis function 



sets |T^, the method of pole (Mittag-Leffier) expansion ||T6|,|T5l of wave functions. Green 



function, scattering amplitudes, etc. 

We feel, that it is very natural to use the so-called J-matrix method in the studies of 
the continuum spectrum effects in nuclear structure or in the low-energy nuclear scattering 
calculations. The J-matrix method was initially proposed for atomic problems fl^jl^ and 
shown to be one of the most efficient and precise methods in calculations of photoionization 



19|,|20[| and electron scattering by atoms |21]. In nuclear physics the same approach has 



been developed independently p2|-p4| as the method of harmonic oscillator representation 
of scattering equations (HORSE). The HORSE method has been successfully used in various 
nuclear applications, e.g., nucleus-nucleus scattering has been studied in the algebraic version 
of RGM based on HORSE (see, e.g, the review papers p5|j26[1 ); the effect of A and neutron 



decay channels in hypernuclei production reactions has been investigated in refs. P7|,p5 



the generalization of the HORSE method ||29| to the case of few-body channels within the 



democratic decay approximation has been used in the cluster model calculations of monopole 
excitations in ^^C pO] and of the halo and soft dipole mode properties in ^^Li [R^ and ^He 



(see also review papers [^,|^), in the study of double-A hypernuclei [^, etc. 
The HORSE method is very attractive because the continuum spectrum wave function 
is represented as a sum of oscillator functions, i.e., in the same manner as in conventional 
variational nuclear structure calculations, say, in shell model or cluster model RGM calcu- 
lations. The wave functions, S'-matrix and other scattering characteristics are obtained by 
algebraic methods that are very efficient in calculations. As a result, it is possible with the 
use of the HORSE method to treat problems with a large number of channels with moderate 
computer facilities (e.g., in J- matrix atomic calculations |T9|j20|l up to ~ 70 channels have 
been allowed for). With the HORSE method it is possible to perform the calculation in the 



complex energy plane and to locate S'-matrix poles numerically [^,^,^,^,0. The loca- 



tion of S'-matrix poles is important not only for calculation of resonance properties p0| , |36 
but also for bound states, e.g., for atomic systems [^ or for weakly-bound nuclear states 
|3T|,|52| the location of the S'-matrix poles improves significantly the variational results for 



binding energies. The expansion of the continuum spectrum wave function in oscillator 
function series involves, of course, an infinite number of terms. However, usually only the 
first few terms of the expansion are important in calculation of matrix elements of opera- 
tors of physical observables. Nevertheless, Nhuj terms with large values of A^ can be easily 
calculated in the HORSE method, and in some special cases these terms become important. 



e.g., in calculation [^ of electromagnetic transition probabilities B{EX;g.s. -^ continuum) 
for weakly bound systems terms up to 4000^a; have been allowed for. 

Below we shall discuss correlations between the HORSE method and the well-known 
i?-matrix and P-matrix approaches (the formal equivalence between the J-matrix method 
and the Feshbach method [Q has been proved in ref. [Q). We shall derive exact expressions 
for R- and P-matrices within the HORSE formalism. These expressions can be used for 
calculation of R- or P-matrices within any variational approach based on the harmonic 



oscillator expansion of the wave functions, e.g., in shell-model calculations, and thereafter 
the well-developed technique of i?-matrix and P-matrix theory (see, e.g., [0,0) can be used 
for calculation of scattering phase shifts, resonance energies and widths, etc., within the 
variational approach. We shall show that within the HORSE formalism a discrete analogue 
of the P-matrix can be introduced that is equivalent to the P-matrix in the quasiclassical 
limit. 

The construction of the P-matrix within the HORSE formalism gives the possibility 
to formulate a very efficient method for introducing the Coulomb asymptotics within the 
HORSE theory. The only method for treating Coulomb interaction within the HORSE 
approach has been developed by the Kiev group |^. The method of ref. [^ involves 



diagonalization of large matrices even in a single-channel case that makes it unusable for 
application to problems with a large number of channels. Our method is free from these 
shortcomings and can be used for high-accuracy calculations of scattering characteristics of 
multichannel systems of charged particles. 

In this paper we discuss binary channels only. 

The paper is organized as follows. Basic equations of the HORSE approach are presented 
in Section 2. P-matrix and P-matrix are derived within the HORSE formalism in Section 3. 
In Section 4, we propose the method of treating Coulomb interaction within the HORSE 
formalism and demonstrate the accuracy of the proposed technique. Section 5 is devoted to 
the generalization of the above results to the multichannel case. Summary is presented in 
Section 6. 

II. POTENTIAL SCATTERING IN THE HORSE METHOD. 

We sketch in this section the basic equations of the HORSE method. We restrict the dis- 
cussion in this section to the simplest case of potential single-channel scattering of uncharged 
particles. 

We use the conventional partial wave expansion of the wave function. 



A'jr 

*k(r) = ^j:t^e^^''ui{k,r)Y;^{n^)YU^^), (1) 

l,m 

where k is momentum, ^2^ and Q^ are angular variables in momentum and coordinate spaces, 
respectively, Yim{^) are the usual spherical functions, 6i is the phase shift in the partial wave 
labeled by the angular momentum /, star is used to denote the complex conjugation. The 
partial amplitudes, ui{k,r), are the eigenfunctions of the radial Schrodinger equation, 

H'ui{k,r) = Eui{k,r) . (2) 

We normalize the partial amplitudes ui{k, r) in such a way that the flux associated with the 
wave function (|l]) is equal to unity. 

Within the HORSE formalism, the partial amplitudes ui{k,r) are expanded in infinite 
series of the harmonic oscillator functions. 



ui{k,r) = J2 ani{k)Rni{r) ■ (3) 



n=0 



Ineq. (|), 



^•'"■' - (-^'" SS^Sy (^)' ^^^ (-fi) ^"" © • '^' 



L"(x) is Laguerre polynomial; the oscillator radius tq = Jh/ jjiu is the only parameter of the 
oscillator basis @; fiu is the spacing between the oscillator levels; /x is the reduced mass. 

The functions ani{k) entering eq. @ represent the wave function in the oscillator repre- 
sentation. They obey the infinite set of algebraic equations 

oo 

J2{Hl^,-5nn'E)an'iik) = 0, (5) 

n'=0 

where the matrix elements of the Hamiltonian in the oscillator basis H^^^, = T^„, + V^^,, and 
T^„, and V^„/ are the matrix elements of kinetic and potential energy operators, respectively. 
The matrix elements of short-range nuclear potentials V^^, — i> in the limit n and/or 
n' —>■ oo, while the kinetic energy is represented by a tridiagonal matrix, 

^nn' = a \n-n' \>1 , 

TL = ^{2n + l + 3/2) , (6) 



TUi,n = Tln+1 = - ^\/(^ + l){n + l + 3/2) 



with non-zero matrix elements T^„ and T^ „_|_^ increasing linearly with n for large values of 
n. Thus, the potential energy matrix can be truncated, and actually we shall use in what 
follows instead of V the potential energy V with matrix elements 

Vl, if n and n' < N; 

VL' = { ~ (7) 

if n and/or n' > N . 

The truncation (^ is the only approximation of the HORSE method. Note, that the kinetic 
energy matrix is not truncated within the HORSE approach in contrast to the conven- 
tional oscillator-basis variational methods which involve the diagonalization of the trun- 
cated Hamiltonian matrix -ff^„/. Thus, the HORSE method is equivalent to the method of 
high-rank separable approximation of the potential. 

We shall introduce in the oscillator representation the interaction region spanned by the 
functions @) with n < N and the asymptotic region spanned by the functions @) with 
n > N. In the asymptotic region, the wave function in oscillator representation, a'^f{k), 
obeys the three-term recurrence relation 

TL, „_i alU, lik) + (TL - E) aim + T^ „+i a^+i, m = , (8) 

as is easily seen from eqs. (^), (|) and (|^). Eq. (H) has two linearly independent solutions, 
Sniik) and Cni{k), and a'^iik) can be expressed as a linear combination of these solutions: 

aZ{k)=cos5iSni{k)+siYi5iCni{k), n>N. (9) 

It is convenient to use the following set of linearly independent solutions Sni{k) and Cni{k) 



US 



S.m ^ /=^^(^^o)- exp (-^) L-nerl) , (10) 

Cum = ^^^^ ,/ "^°"' ,(kr,)-^ exp (-^^] 

"'^ ^ r(-/ + i/2) V^r(n + / + 3/2) ^ °^ ^\ 2 J 

x$(-n-/-l/2, -/ + 1/2; A;V^) , (11) 



where the velocity v = hk/jj,, and <l>(a, b; z) is a confluent hypergeometric function [HD 
The solutions (0) and (pUD are defined in such a way that 



Y,SUk)Rni{r) 



n=0 



3i{kr) 



(12) 



and in the limit r — >^ cxo 



°° k 
Y, Cni{k)Rni{r) -^ y= ni{kr) . 

n=0 V^ 



(13) 



Here j/(x) and ni(x) are spherical Bessel and Neumann functions [^,0, respectively. Thus, 
the adopted normalization of the wave function (|l|) is assured, and 6i entering eq. (|]) appears 
to be just the phase shift in the partial wave corresponding to the angular momentum /. 

The functions (|T0|) and ( pJ]) can be easily calculated. The function Sni{k) is seen to be 
the harmonic oscillator function in the momentum representation. 

We shall need in what follows asymptotics of the functions (p!0| ) and ( [Tl| ) in the limit 
n ^ oo: 



Sni{k) ^ 2kroJ- (n + 1/2 + 3/4)3 ji{2kro^n + 1/2 + 3/A) 



and 



— {n + 1/2 + 3/4)~3 sin[2A;ro ^n + 1/2 + 3/4 - 7:1/2] 



(14a) 
(14b) 



Cni{k) ^ -2krJ— {n + 1/2 + 3/4)^ ni(2krQ^n + l/2 + 3/A) 



{n + 1/2 + 3/4)^3 cos[2A;ro ^n + 1/2 + 3/4 - 7x1/2] , 



(15a) 
(15b) 



that can be easily obtained from the asymptotics of the confluent hypergeometric functions 
| 40| and Laguerre polynomials |41 . 

The so-called Casoratian determinant, /C^(C, 5) = Cn+i,i{k) Sni{k) — Cni{k) Sn+i,i{k), 
plays the same role in the theory of second-order finite difference equations like eq. (|^) (see, 
e.g., [43, 44 1) as Wronskian in the theory of second-order differential equations. It is easy to 
show that A = T/j ,„_,_]^/CJj(C, 5*) is n-independent. Thus, A can be easily calculated in the 
limit n — > cx) using asymptotics ( |14b|) and ( |15b|) . As a result, we get: 



T: 



n,n+l 



/C^(C, 5) 



h 

2 



(16) 



In the interaction region n < N, the set of algebraic equations (j^) due to (H) and (^ 
takes the following form: 

Y: (i^L' - Snn'E) a^n{k) = -<5„^ T^^^^^ a^^i,^/^); n = 0, 1, . . . , iV . (17) 

n'=0 

The functions a^iik) satisfying eqs. (^) in the interaction region n < N can be expressed 
in terms of a'^_^i i{k), 

ani{k) = GnN a^^+i^iik) , (18) 

and the matrix elements Qnn' can be calculated by the formula 

IXn 7A?i' rpl 



^-' = -E^^^^^n;n'+l. (19) 

A=0 



Ex — E 



Here Ex and jxn are eigenvalues and corresponding eigenvectors of the truncated Hamilto- 



nian matrix -ff^„/, n,n' = 0,1, ... ,N. 



The function a^iik) should fit both eqs. (H) and (|T^. Thus, we can express both a^iik) 
and aJv+i i{k) in terms of the phase shift 5i using eq. (^), and insert these expressions in eq. 
(^). As a result, we obtain a simple formula |TE[ for the phase shift 5/: 



■- _ SNljk) - QNNSN+l,l{k) , . 

So, within the HORSE method, one should first diagonalize the truncated Hamiltonian 
matrix if^„;, n,n' = 0, 1, . . . , A^, i.e. one should find eigenvalues Ex and eigenvectors •yxn- 



Next, the matrix Qnn' and the phase shift 6i are calculated using eqs. (0) and ( pO]) while 



Sni{k) and Cni{k) entering eq. (pO]) are calculated using (|TOD and ([TT]). The functions aniik) 
defining the wave function in the oscillator representation, are to be obtained by eq. (P) in 
the asymptotic region n > N and by eq. (IT^) in the interaction region n < N. 



All the above expressions can be generalized to the case of multichannel scattering |T^ 
or to the case of true few-body scattering (democratic decay approximation) [^. The 
convergence of all of the physical observables in nuclear calculations is usually obtained if 
the interaction region for each channel is spanned by 5-10 oscillator functions (see, e.g.. 
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The truncated Hamiltonian matrix Hnn' is actually the same matrix that is used in con- 
ventional variational approaches with oscillator basis, e.g., in the shell-model calculations. 
The HORSE method is a natural generalization of these variational approaches on the case 
of continuum spectrum states. The calculation of the matrix elements of the truncated 
Hamiltonian matrix Hnn' and its diagonalization is the most complicated part of the cal- 
culations within the HORSE approach. An important point concerning the applications of 
the HORSE method is that the diagonalization of Hnn' should be done only once, thereafter 
the negative eigenvalues E\ can be interpreted as the energies of the bound states and the 
corresponding eigenvectors 7a„ can be used instead of a„; in the expansion (H) for the con- 
struction of the bound state wave functions, while continuum spectrum wave functions in a 
wide range of positive values of energy E are easily calculated by means of simple formulas 
(I), (0), (0), (0), (0) and (H). 

III. DISCRETE ANALOGUE OF THE P-MATRIX AND NATURAL CHANNEL 

RADIUS. 

The P-matrix is identical to the inverse i?-matrix and is defined 0J^ as the logarithmic 
derivative of the wave function at the channel radius r = b: 



P = R 



_,_HMk,r] 



(21) 



ui{k,b) 

The partial amplitudes ui{k, r) = rui{k, r) are often used in applications. We can define also 
P-matrix P by the expression (^) with Ui{k, r) instead of Ui{k, r). It is easy to see that 

P = P - 1 . (22) 

The channel radius 6 is a parameter of the P- and P-matrix theories; b divides the whole 
coordinate space into two parts: the "interaction region" r < b and the "asymptotic region" 
r > b, i.e. the channel radius b plays the same role as the truncation boundary A^ in the 
oscillator space of the HORSE method. In the asymptotic region r > b, the wave function 
is of the form: 

9 



ui{k,r) = COS Si Fi{k,r) + sin 6iGi{k,r) . (23) 

If some long-range interaction, say, Coulomb interaction, cannot be neglected in the 
asymptotic region r > b, than Fi{k, r) and Gi{k^ r) are respectively the regular and irregular 
solutions of the radial Schrodinger equation involving the corresponding potential. In the 
simplest case, when the interaction between particles can be neglected in the asymptotic 
region r > b, Fi{k,r) = -^ji(kr) and Gi{k,r) = — y=ni{kr). If, additionally, the channel 



radius b is large enough and spherical Bessel and Neumann functions can be replaced by 
their asymptotics at r > 6, then in the asymptotic region 

ui{k, r) = — - sm( kr - — + Si) . (24) 

ry/v ^ 2 ^ 

In this case, from eq. (|2l|) it follows, that 

P = R-^ = kb cot^kb - — + Si) , (25) 

while in the general case, we have 

p _ ^-1 _^ FI{k,r)+ tan SiG'iik,r) 
Fi{k,b)+teinSiGi{k,b) ' 



d 



r=b 



where F/(fc, b) ^ -^Fi{k, r) and G'lik, b) ^ -^Gi{k, 

r=b 

According to the general theory of ref. [0, in the interaction region r < b the wave 
function ui{k, r) is expanded in a complete set of functions {wnz(^)} matching some boundary 
condition at r = 6, e.g., in the P-matrix approach f„;(r) satisfy the boundary condition 
Vni{b) = 0. Using the formalism described in detail in refs. 0,^, P- or i?-matrix can be 
calculated solving the Schrodinger equation in the functional space spanned by the functions 
{fJni{r)}, and thereafter the phase shift Si can be found using eq. (^) or (|26D. The well- 
developed formalism of R- or P-matrix theory can be also used for calculations of resonance 
positions and widths. 

It is seen that the HORSE formalism has much in common with the P-matrix and P- 
matrix ones. All the above mentioned formalisms involve the division of the whole space 

10 



into the interaction and asymptotic regions. Nevertheless, within the HORSE approach one 
divides the functional space, while within the R- and P-matrix approaches the coordinate 
space is divided into two regions. The HORSE formalism involves the oscillator basis that 
differs significantly from the basises used in the i?-matrix or in the P-matrix theory, and 
there is no formal equivalence between the HORSE method and the P-matrix method or 
the P-matrix method. The HORSE formalism is very attractive because oscillator basis is 
natural in nuclear structure studies, while P- and P-matrix formalisms involve basises that 
are more artificial for nuclear physics applications. 

Applications of the P-matrix formalism |]3|j45[| refer to a very important property of the 
P-matrix: the P-matrix has poles at the energies that can be identified with the energies 
of the so-called 'primitives' [Q], i.e., of the eigenstates obtained with an artificial boundary 
condition f„/(&) = 0. This feature correlates scattering phase shifts, resonance positions and 
widths with the structure of the spectrum of the Hamiltonian in the interaction region in 
the P-matrix approach. The boundary condition of the type Vni{b) = is good for quark 
bag models 0,0, but it is dubious for nuclear structure applications. If one starts from 
the nuclear shell model, than it is not clear what value should be attributed to the channel 
radius b and what is the correlation between the primitives and the shell model states. Thus, 
it is interesting to derive the P-matrix within the HORSE formalism. 

Let us suppose that we can neglect the interaction at distances r > b. Then for calculation 
of the P- or P-matrix one can use eq. ( p6D or eq. (^). Substituting tan 5/ in eq. ( ]26|) by 



the expression (pO]), we obtain the following expression for P- and P- matrices within the 



HORSE formalism in the general case: 

^ ^^-1 ^ Cr^iimikb) + Sr,i{kX{kb) - g^r, [C^-,iAm{kb) + S^+,j{kX{kb)] ^27) 
Cm{k)ii{kb) + Sm{k)ni{kb) - Gnn [CN+i,i{k)ji{kb) + SN+i,i{k)ni{kb)] 



d 



d 



where j'likb) = -^ji{kr) ^ andn'i{kb) = -^ni{kr 



r=b 



In the case when the approximation 



r=b 



can be used, expression (pTf ) can be simplified and it takes the following form: 
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Cm (k) + tan (fcfe - f ) Sm (k) - Gnn C^+i, / ( fc) + tan ( fc6 - f ) S^+i, i (k) 

— rvU 



Cm{k) tan (fc6 - f ) - Sm{k) - Gnn 



CN+i,i{k) tan [kb - f ) - 5^+1, K^) 



1. (28) 



Eqs. ( pTl) and (^) are very general and can be used within the HORSE formahsni for any 
value of the channel radius b large enough to neglect the potential energy at r > 6. These 
equations can be used for calculation of P- or i?-matrix within the shell model approach 
with the subsequent use of the well-developed R- or P-matrix formalism for calculations of 
phase shifts, properties of resonances, etc. 

It is seen, that the P-matrix poles do not coincide with the eigenenergies of the truncated 
Hamiltonian matrix that are actually the poles of Qnn as evident from the definition (P^. 



Thus, generally the oscillator shell model states cannot be associated with the poles of the 
P-matrix. To put the P-matrix poles in the one-to-one correspondence with the shell model 
states, we should remove Qnn from the denominator of eq. (|27D . This can be done if the 
channel radius b fits the following equation: 

jijkb) ^ _ SN+i,i{k) 

ni{kb) Cn+iM ■ ^ ' 

Eq. (p9D has an infinite number of energy- dependent solutions 6j. Nevertheless, in the 
quasiclassical limit A^ -^ oo, one of the solutions becomes energy-independent. This solution, 
6o, is easily found substituting SN+iAik) and CN+i,i{k) in eq. ( p9D by their asymptotics ( |14aj ) 
and ( |15a| ): 



&o = 2royiV+l + I. (30) 

Eq. (^) may be reduced to a much simpler form ii ji{kb) and ni{kb) are replaced by their 
asymptotics that is equivalent to the use of the approximation (|2^) for the wave function in 
the asymptotic region and eq. (|28|) for P-matrix: 

This equation can be easily solved, and we obtain: 

12 



7tI 1 /•S'at.li i(k)\ in , . 

^- = ^+T:^^"^^^ / M +17' ^ = 0,±l,±2,... . 32 

2k k \CN+i,i{k)J k 

The solution ( |5(]| ) corresponds to i = that is easily seen substituting SN+i,i{k) and 

CN+i,i{k) in eq. ( |5^ ) by their asymptotics ( |14b| ) and ( |15b| ). It is seen that all the rest 



solutions bi -^ bo in the limit k ^ oo. 

We shall refer to the energy-independent channel radius defined according to eq. (^) as a 
natural channel radius bo. The natural channel radius 6o coincides with the classical turning 



point, r^+i = 2ro JN + 1/2 + 7/4, corresponding to the oscillator function R]\[j^ij{r), i.e. 
to the first function outside the interaction region in the oscillator representation used for 
the construction of the truncated Hamiltonian matrix Hnn'- 

It is worthwhile to note, that though eq. ( |30D has been obtained using the asymptotics 



|14a| ) and ( |15a| ), it fits well one of the solutions of eq. ( pQ] ) even for very small values of 



b — b 
N. This is illustrated by fig. 1 where the relative deviation, — ; — , of the natural channel 



radius 6o in various partial waves from the exact solution b of eq. (EUI) is plotted vs energy 



E. The calculations are performed with hu = 18 MeV and reduced mass corresponding to 

b — b 
the neutron scattering by A = 15 nucleus. The absolute value of — ; — is seen from the 



figure to increase linearly with energy. However, even in the limiting case A^ = (only one 

b — b 
oscillator function present in the interaction region), the relative deviation — - — does not 

exceed 8% only. 

Using eqs. (|l^) and (pT]) we obtain for the P-matrix corresponding to the channel radius 

bi matching eq. (^) the following expression: 

P= r ni{kbi)CN+i,i{k) 

n 

X {Ukh)Cm{k) + Sm{k) - Gnn [Ukh) C^+iAk) + ^TV+i,^^)]} , (33) 



where ji{kb) = i[{kb) / n'i{kb) and ni{kb) = n'i{kb) / ni{kb) . In the approximation (plD , eq. (^ 

reduces to 

2kb r' 

P = '-^'^^[CN+iMCNi{k) + SN+i,i{k)Sm{k) - QNN[cl+i,i{k) + 5^+i,,(A;)]}-l. 



(34) 
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Further simplifications of eq. ( ^41) can be introduced using the asymptotics (|14b| ) and 
( |15b| ) of the functions Cni{k) and Sni{k). After some straightforward algebra, in the most 
interesting case b = Bq we derive a very simple expression connecting P- matrix and Qnn'- 



P = 2 J(iV + l) (iV + / + |) [(3-Gnn) - 1 , (35) 



where 



^ 2N + 1 + 1/2Y 
^=i 2iV + / + 3/2 J ^°^ 



2^roh/iV+i + I-JiV+i + ^ 



(36) 



Obviously, /5 -> 1 in the limit of large N, and ^(A^ + l)(iV + / + 3/2) = N + 1/2 + 5/4 + 
0{N"^). Thus, we can rewrite eq. (PH) as 



P = 2{N + l + -^[l-gM^)-l, (37) 

or as 

2Ar + / + ^) ^^^±Mi^l^^ - 1 (38) 

using eq. (|TB]). 

Equation (133) obtained in the quasiclassical limit N ^ 00 can serve, however, as a good 
approximation for the exact P-matrix even for very small values of N. As is seen from fig. 2a, 
even in the case A^ = 1 the P-matrix calculated using eq. (|38|) is very close to the exact one 



calculated by eq. (pT]). However, there is a small discrepancy between the positions of the 
poles of the P- matrices and, as a result, between their values in the vicinity of the poles. 
As the number of oscillator functions in the interaction region increases, the discrepancy 
disappears and in the case A^ = 9 the low-energy poles coincide (see fig. 2b). However, there 
is still a small discrepancy between the poles at higher energies that arises from the small 
increase with energy of the channel radius b fitting eq. (^) while the natural channel radius 
60 is energy-independent (see fig. 1). 

It is interesting to compare eq. (^) with eqs. (^T])-(^). The first term in eq. ( P5[ ) 
looks like a finite-difference analogue of the P-matrix P with the wave function in oscillator 
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representation aN^ij{k) and finite difference Aqn+i, i{k) = a^^ij^k) — aNi{k) playing the 



roles of the wave function ui{k,b) and its derivative ^ui{k,r) , respectively. 



r=b 



To make this analogy more transparent, we note that the range of an oscillator function 
Rni{r) is characterized by its classical turning point, r'^ = 2rQJn + 1/2 + 3/4. Thus, in the 



definition of the discrete analogue of the P- matrix we should use instead of ■^ui{k, r 
finite-difference analogue, 

AciN+ijik) _ aN+ij{k) - am{k) 



r=b 



its 



A<+i 



^Af+l ^Af 



aN+i,iik) - am{k) 



ro 



Ur ^ 5 r 
2 4 



ON 



-1 



(39) 



Contrary to the coordinate space, we have two characteristic lengths in the harmonic os- 
cillator space instead of a single length parameter b: the range of the interaction space, 
rff, equal to the range of the last basis function R^iij-) included in the interaction re- 
gion, and the range of the first basis function Rjq^iiir) included in the asymptotic region, 
rffj^i- Defining the discrete analogue of the P-matrix, Vn-, we replace 6 by a symmetri- 
cal combination of both lengths. 



^N ^N+1 



2ro[{N + 1/2 + 3/4) (AT + 1/2 + 7/4)]!/^ 



2roJN + 1/2 + 5/4 [1 + 0{N ^)]. As a result, for the discrete analogue of the P-matrix P 
we obtain the following expression: 

1 AaN+i,i{k) 



^''^'''-"'"■^-^^^^u^)^?^ 



2N + 1 + 



N+l 

5\ aN+i,i{k) - amik) 



(40a) 
(40b) 



27 aN+i,i{k) 

(terms proportional to A^~^ and higher-order corrections are omitted in eq. (|40b| )). So, we 



have shown that the first term in eq. (|38|) is really the discrete analogue Vn of the P-matrix 
P. 

The P-matrix P involves the partial amplitudes Ui{k, r) while the P-matrix P involves the 
partial amplitudes ui{k, r) = -ui{k, r). By analogy, in the harmonic oscillator representation 
we define functions ani{k) = -^aniik). Replacing aniik) in eq. ( |55D by r'^ani{k) and omitting 



terms proportional to A^ and higher-order corrections, we rewrite eq. 



as 
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It is easy to show that the right-hand side of eq. (0) is just the discrete analogue of the 
P-matrix P defined as 

aN+i,i{k) Ar^^i 

Concluding this section, we have derived expressions that can be used for calculation of 
R- or P-matrix in any variational approach based on the harmonic oscillator expansions, 
e.g., in the shell model approach. We have shown that if only the channel radius is equal to 
the natural channel radius or any other radius 6j fitting eq. (129|) or eq. (^2t), the variational 
energies coincide with the P-matrix poles; in the quasiclassical limit, the P-matrices P and 
P for the natural channel radius reduce to their discrete analogues Vn and Vn, respectively. 
These results link the inner structure of the system with scattering characteristics and do 
not depend on the interaction in the system, on the structure of the energy spectrum, etc. 

IV. COULOMB INTERACTION WITHIN HORSE FORMALISM. 

A. Theory 

We have discussed above only scattering on short-range potentials and have not allowed 
for the Coulomb interaction. Till now the only method of treating the scattering of charged 
particles within the HORSE approach is the one that has been proposed by Kiev group in 
ref. [39|. Matrix elements of the Coulomb potential in the oscillator representation, V^°i^\ 



decrease very slowly as n or/and n' increases. Thus, following the lines of ref. [^, to 
allow for the Coulomb interaction within the HORSE formalism one should increase the 
truncation boundary A^ up to the values A^ ~ 50. In the case of long-range Coulomb 
interaction, analytic expressions like (10) and ([TT|) for the linearly independent solutions of 



the Schrodinger equation in the asymptotic region n > N are unknown, and to calculate 



the solutions one should use asymptotic expressions [^. Using the method of ref. [^, it 
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is possible to calculate with high accuracy scattering observables in the case of few-channel 
scattering of charged particles, but application of the method to the case of a large number 
of channels is not feasible. 

The P-matrix discussed in the previous section can be used for formulation of a ver- 
sion of the HORSE formalism allowing for the Coulomb interaction that is free from the 



shortcomings of the method proposed in ref. ||3^ 



Let us consider the scattering in the system with the interaction described by the poten- 
tial V = v^^«c' + V"^°''', where 1/^"^' is a short-range nuclear potential, and y^°"' = ZiZ2e^/r 
is the Coulomb potential. In this case instead of (|l]) the wave function is expressed as 

^k(r) = -j:^^e^'^'''~''^Mk,r)Y;^{Q^)YU^^) , (43) 

l,m 

where t]i = argr(l + I + i() is the Coulomb phase shift, and Sommerfeld parameter ( = 
ZiZ2e^lJi/k. We introduce the channel radius h that is supposed to be large enough to neglect 
the nuclear potential y^"^' at the distances r > b, i.e., we suppose that b > RnucI where 
RnucI is the range of y^"^'. In the asymptotic region r >b, the partial amplitudes ui{k,r) 
are expressed by eq. @ with Fi{k,r) = ^ ^Fi{C,kr) and Gi{k,r) = - ^ ^Gi{C,kr), 
where regular Fi{(,kr) and irregular Gi{(,kr) Coulomb functions are defined according to 
ref. [0. The P-matrix, P^ , accounting for the Coulomb asymptotics, satisfies the general 
expression (^61) . 

Let us introduce an auxiliary short-range potential V^^ by cutting Coulomb potential 
yCoui g^^ ^Yie point r = b, i.e. 

{Y = Y^"^^^ -\- yCoMi r <■ b 
' " ; b>RNuci. (44) 

0, r>b 

The continuum spectrum wave function corresponding to the potential V^^ is given by (|l]) 
with 6i = 6f^, where 6f^ are the /-wave scattering phase shifts for the potential Y^'^. In the 
asymptotic region r > b, the corresponding partial amplitudes Ui^{k,r) are of the form: 

uf\k,r) = cos5fVK^,0 + sm6i^^gi{k,r) , (45) 
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where fi{k, r) = -y=ji(kr) and gi{k, r) = — ^ ni{kr). The Schrodinger equation for the short- 
range potential V^^ can be solved within the HORSE formalism, and the corresponding 
P-matrix, P^'^, can be calculated by eq. (^Tl) . 

An important point is that the P-matrix, i.e., the logarithmic derivative of the wave 
function ui{k,r) at the point r = 6, is dictated by the interaction region r < b only, and 
does not depend on the potential in the asymptotic region r > b. Thus, P-matrices P*" and 
pSh corresponding to the potentials V and V^^, respectively, should be equal at the point 
r = 6, 



P^ = P 



Sh 



(46) 



Substituting P^ and P^^ in eq. p6|) by the right-hand sides of eqs. ( ]26|) and (pTf), respec- 
tively, we derive the following expression for the phase shift 5i corresponding to the potential 

T/ '\rNucl _j_ '\rCoul. 



tan 5/ 



Cm{k) - g'^%CM+iAk) WbUuFi) + 



Cmik) - Gm^CN+iAk) Wbiji,G, 



SNi{k) — GN]\rS]\f+i^i{k) 



SNiik) 



g^%SN+iAk) 



Wt,{ni,Fi] 



Wbini,Gi) 



(47) 



Here, g^% is defined by eq. (|T^) and corresponds to the Hamiltonian with the auxiliary 
short-range potential V^'^'*, the quasi- Wronskian Wb{ji,Fi) = 

{£[Ukr)]Fi{C,kr) - Ji{kr)£Fi{C,kr)}l^^, and Wk{ni,Fi), W,Ui,Gi) and W,{ni,Gi) are 
expressed similarly. 

Eq. (|47|) defines within the HORSE formalism the phase shift 6i in the case of scattering of 
charged particles. It involves diagonalization of the Hamiltonian with short-range interaction 
ySh QT^ij Thus, moderate values of the truncation boundary A^ can be used; as a result, 
the corresponding algorithm appears to be much more effective as compared to the method 
of ref. p|. 

Taking into account eq. (pOD, it is easy to obtain from ( ^7D the expression for the phase 
shift Si in terms of the phase shift S^^ corresponding to the auxiliary short-range potential 

ySh. 



Sh 



, e W,{ji,Fi) -W,{ni,Fi) tan 6 



W,{ji,Gi) - Wi,{ni,Gi) tan 6, 



Sh 
I 



This expression has been proposed in ref. [|7] for momentum-space calculations of charged 
particle scattering. 

The partial amplitude uf^{k,r) corresponding to the potential V^'^ obtained by the 
HORSE approach, perfectly matches in the interaction region r < b the required wave 
function ui{k, r) corresponding to the potential V = y^^f^^j^yCoui^ |^^^ uf^{k, r) and u/(/c, r) 
differ at larger distances where the asymptotics ( p3| ) with Fi{k,r) = ^ J^Fi{(,kr) and 
Gi{k,r) = — ^ vl;^KC) ^^) should be used instead of (^51). However, typically only the 
interaction region is of importance for calculations of electromagnetic transition probabilities 
and other observables. At the same time, it should be taken into account, that asymptotics 
are only significant in calculation of normalization factors for continuum spectrum wave 
functions. Thus, before calculating matrix elements of operators of physical observables. 



u 



Sh 



{k, r) should be renormalized. The renormalization factor A/" can be easily derived from 



the equation ui{k,b) = N'uf^{k^h) where Ui{k,r) and Ui{k,r) are given by (^) and ( p3| 
respectively: 



cos6iFi{C,kb) - sin 6iGi{C,kb) 
cos 5i^ ji{kb) — sin 5i^ni{kb) 



■^^ ~ „^„ ^Sh ■ (Uh\ ^;„ xSh^ n^l.\ ■ y^^) 



B. Numerical illustration. 

We illustrate the applicability if our method for treating Coulomb interaction within the 
HORSE formalism by calculations of the single-channel proton-nucleus potential scattering. 
We use Woods-Saxon potential with conventional surface spin-orbit term as the nuclear 
potential, V^^'^K Parameters of \/^"^' and of the Coulomb potential V'-'"^'' taken from ref. [^ 
correspond to the p-^^N scattering. Most of the result presented below were obtained with 
the oscillator level spacing parameter hu = 18 MeV. We made use of Lanczos smoothing 



of potential energy matrix elements [^ that improves the convergence. Note that the 
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convergence of the results of HORSE calculations without Coulomb interaction is usually 
achieved at A^ ^ 5 ^ 6. 

Channel radius b used for construction of the auxiliary potential V'^^, is a free parameter 
of the method. It should be taken larger than the range of the Woods-Saxon potential 
RnucI- On the other hand, the truncated Hamiltonian matrix -ff^„/ should carry information 
about the jump of potential V^^ at the point r = b. Thus, b should be chosen less than the 



classical turning point, r^ = 2ro JN + 1/2 + 3/4, of the oscillator function RNi{r), i.e. of 
the function with the largest range in the set of oscillator functions {Rni{r), n < N} used for 
the construction of the truncated Hamiltonian matrix H^^,. Figure 3 gives the s-wave phase 
shifts 6q corresponding to the proton energies E = 2 and 10 MeV calculated with different 
values of b. In our case, RnucI ~ 5.5 fm while hu = 18 MeV corresponds tq ^ 1.5 fm and 
hence N = 9 corresponds r^ ?^ 9 fm. It is seen from fig. 3 that the plot of So{b) has a plateau 
between RnucI and r^ that reproduces well the exact values of 6o. Moderate variations of 
huj cause only very slight changes of the phase shift while the increase of N results only in 
extension of the plateau to larger values of b. The phase shifts 5/ (6) with / > also have 
the plateaus that reproduce the exact values even better than in the case / = 0. 

So, one should assign to the channel radius b a value from the interval RnucI < b < rfj. 
Results presented below were obtained with 6 = 7 fm. 

Plots of the s-wave phase shift 6^ versus proton energy E are presented in fig. ^. The 
exact values of ^o were obtained by numerical integration of the Schrodinger equation in 
the differential form. It is seen that the results obtained in our approach exactly reproduce 
the phase shift np to E = 30 MeV and only the phase shift obtained with A^ as small 
as 4 deviate a little from the exact ones. We note that there is a sharp resonance in the 
s-wave scattering. The resonance is generated by the Coulomb barrier and disappears if 
the Coulomb interaction is turned off. The phase shift in the vicinity of the resonance 
depends crucially on the resonant energy. Our calculations reproduce well the position of 
the resonance. It is interesting that the best description of the resonant energy is obtained 
in calculations with A^ = 4, and hence the A^ = 4 phase shift is closer to the exact one in the 
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vicinity of the resonance than the phase shifts obtained with larger values of A^. We note, 
however, that the height of the jump of the phase shift in the vicinity of the resonance is 
reproduced in calculations with larger values of A^ better than in calculations with A^ = 4. 



For comparison, the phase shifts calculated by the method suggested by Kiev group p9 



are also plotted in fig. B. The thin solid curve with diamonds was obtained with the trunca- 



tion boundaries recommended in ref. |3^: the Hamiltonian matrix was truncated at A^ = 70 
and hence 71 x 71 matrix was diagonalized in this single- channel problem, however the ma- 
trix of the short-range nuclear potential was truncated at A^"^'' = 50. It is seen that the 
Kiev method describes the phase shifts with the same accuracy that ours (note, however, 
small deviation of the Kiev phase shift from the exact one in the low-energy region) but is 
much less efficient in applications because it involves diagonalization of a large matrix. If 
the truncation boundaries of the Kiev method are reduced essentially, the description of the 
phase shift becomes worse. This is illustrated in fig. ^ by dots which present the phase shift 
obtained with A^'^'^ = 8 and A^ = 70, i.e. the short-range potential was truncated to 9 x 9 
matrix while the Coulomb interaction was truncated to 71 x 71 matrix, and therefore in 
calculation of the phase shift it was needed to diagonalize the 71 x 71 Hamiltonian matrix. 
Note that contrary to the Kiev method, we obtain an excellent description of the phase shift 
with A^ = 8, i.e. diagonalizing the 9x9 Hamiltonian matrix. If the matrix of the same rank 
is diagonalized in the Kiev method, i.e. the truncation boundary A^ is reduced by an order 
of magnitude, the method fails in the description of the Coulomb asymptotics of the wave 
functions, of the phase shifts, etc. 

The scattering phase shift is a very important but not the only characteristic of the con- 
tinuum spectrum states. In calculation of electromagnetic transition probabilities and other 
observables, one needs matrix elements of various operators. Within the HORSE formalism, 
the matrix elements are expressed through the standard matrix elements in the oscillator 
basis (^ of the operator of interest and the wave function in oscillator representation ani{k) 
entering expansion (^. The plots of a'^i{k) for / = and n = 0, 1, 2 are given in fig. ^. The 
exact values of a^i{k) were calculated as 
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aliiq) 



oo 

Ui{k,r)Rni{r) r'^dr 





(50) 



where the wave function ui{k,r) was obtained by numerical integration of the Schrodinger 
equation in the differential form. It is seen that the wave function in oscillator representation 
ani{k) is reproduced with the same high accuracy in our method as the phase shift. The 
description of ani{k) in the Kiev method is also approximately of the same accuracy as the 
one of the phase shift. The most important though small enough deviation from the exact 
result is seen only in the vicinity of the resonance and is dictated by the resonant energy 
that differs a little in various approximations. Note that the shape of the resonant curves is 
well reproduced for all values of truncation boundary. The best description of the resonant 
energy is obtained accidentally in the N = 4 calculation, and the A^ = 4 curves are closer to 
the exact results in the vicinity of the resonance than the ones obtained with larger values 
of A^. Note, however, that the integral ^ a\i{k) dE in the vicinity of the resonance is better 
reproduced in A^ > 4 calculations. 

The coordinate space radial wave functions ui{kr) are reconstructed within the HORSE 
formalism as 



M 



ui{k,r) =rJ2(^niik) Rrd{r) . (51) 

n=0 

The plots of ui{kr) are given in fig. ^ It is seen that the wave function deviates from the 
exact one if M = iV even in the case when the truncation boundary A^ is large enough 
to reproduce exactly the phase shifts. The deviation is easily eliminated by setting M 
to be large enough and allowing in the expansion (|5T|) for the wave function in oscillator 
representation in the asymptotic region a^/(fc) with n > N that is calculated by a simple 
formula (^. It is interesting to note that in the great majority of nuclear applications like 
calculations of electromagnetic transition rates and of other observables, allowing for afai^) 
with n > N does not effect the results, and hence the deviation of the wave function ui{kr) 
obtained by (^TJ) with M = N from the exact one that seems to be essential, is really of no 
importance. 
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The p-wave and (i-wave phase shifts and corresponding a^iik) are plotted on figs. ^ and 
^. It is seen that the p-wave and rf-wave observables are reproduced in our approach with 
the same accuracy as the s-wave ones. 

Concluding this section, we developed the Coulomb HORSE approach that is efficient 
in calculations of continuum spectrum wave functions and scattering observables within the 
HORSE method in the case of short range + Coulomb interaction. 

V. MULTICHANNEL SCATTERING. 

In this section we generalize the results of previous sections to the multichannel case. 



In the multichannel case, the wave function is of the form Pq , ^^ 

Expression ( ^2|) can be used both for charged or neutral particles, in the latter case the 
Coulomb phase shifts rji^ = 0. It is supposed that the incident plane (or Coulomb-distorted) 
wave in the case of \l/f^ or the outgoing plane (or Coulomb-distorted) wave in the case of \I^|~ 
is present only in the "physical" channel labeled by the multiindex Fj = {z/, fcr^, sr^, err,}, 
while outgoing waves in the case of \l/f, and ingoing waves in the case of \l/^~ are present 
in all "physical" channels F = {z/, fcr, sr, o"r}. Here, sr and ar are the total spin of the 
colliding particles in the channel F and its projection, respectively; kr is the relative motion 
momentum in the channel F; all the rest quantum numbers in the channel F, in particular 
the ones related to the internal state of the colliding particles described by the internal 
wave function ^uirsrixr), are labeled by the index z/. In the right-hand side of eq. ( |52[ ) 
we use multiindex F = {z/, J, /p, sr, crp} to label channels characterized by a definite value 
of the total angular momentum J, Zp is the orbital angular momentum in the channel F. 
The transformation of the states of the channels F with a definite J into the states of the 
"physical" channels F characterized by the vector of the relative motion momentum kr of 
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the incident or outgoing plane (or Coulomb-distorted) wave, is performed with the help of 



(see, e.g., M]) 



X {lrmrST(TT \JM){lrmr^sr,(Jr, \JM) 



(53) 



where {Imsa \JM) stands for a Clebsch-Gordan coefficient, mp is the projection of the 
orbital angular momentum Ir in the channel F, and M is the projection of the total angular 
momentum J. 

The relative motion in the channel F is described by the radial wave function 
pL N (/i;r,rr). In the case when the incident plane wave corresponds to the channel Fj 



u 



and outgoing spherical waves are present in all channels, the radial channel functions are of 
the form in the asymptotic region rr > &r : 



u 



Ji+) 



r{+) 



r{-) 



(54) 



where Srv is a matrix element of the S-matrix, 

Hi^\k,r) =Gi{k,r) ±tFi{k,r), (55) 

and Fi{k,r) and Gi{k,r) are defined in sect. III. In calculations of photodisintegration cross 
sections one needs radial channel functions u-p/p\ {kr,rr) = [ui't-pK ikr,rr)) characterized 
by the outgoing plane (or Coulomb-distorted) wave in the channel Fq and ingoing spherical 
waves in all channels. Below we shall discuss the radial channel functions u-p,-p{ {kr, rp) only; 
of course, the same formalism can be applied to the functions Mpfr ) (^r,''"r)- 

It is convenient to introduce radial channel function matrix (we suppose that the number 
of open channels is M and use square brackets to denote M x M matrices in the channel 
space) , 



(56) 



M^(+) {k, r) 


— 


^iw (^2, ^2) 


^2(2) («^2,r2) ■ 


■ %(M)(^l'^l) 

■H+) n \ 

■ «2{A/)(^2,r2) 




1L,l,l (knj.rAr) 


v:{it\(k...r..r) ■ 


■ V:{it},.(k.r.r.,r) 
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A natural generalization of eq. ([21| ) leads to the following definition of a multichannel P- 
matrix P8|j3[1: 



[P] = [h] u'^+\k,h) u'^+>{k,b) 



,J(+)( 



[h] \ [i7(-) (A;, 6)1 ' - [i^(+) (A;, h)] ' [S]\\ W'^ {k, 6)1 - [i7(+) (A;, 6) 



(57a) 
. (57b) 

The symmetry of the multichannel P-matrix [P] is discussed in the Appendix. In 
eqs. (p7i|)-([57ED, 



«■'(+) (A;, 6) 
/ 



«iaV(^i'^i) 



dri 



«i(2V(^i'^i) 



Ji+) 



^(1)^(^2,^2) 



dro 



ri=bi 



r2=b2 



dri 
ui\2\ (^2, ra) 



"2(2) 



dro 



ri=bi 



r2=b2 



Ji+) 



'"m(i) v^m-i'J'm) 



\ 



dr 



M 



'^M(2) i^M, ^m) 



rAi=bA 



dr 



M 



rM=bA 



[S] is the S'-matrix, the diagonal channel radius matrix 



/ 



[b] 



h 



\ 



y hu J 



and diagonal matrices 



^ Hlf^{hM) 



[^(±) (A;, 6) 



^^^(^2,62) 



and 



Wi{m)(^1'^i) 



dri 
ui(M) (^2, rs) 



''2(M) 



dro 



ri=bi 



r2=b2 



Ji+) 



"^MiM) i^M, ^m) 



dr 



M 



rM=bA 



H\l\kMM) 



(58) 



(59) 



(60) 
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/ rfif w (^^^ ^^) 



H^^^k,b)]' = 



dr^ 



ri=6i 



dHi^^ ih, r,] 



dro 



r2=b2 



dHit^ {kM, Tm) 



dr 



M 



rM=bi\ 



(61) 



Within the multichannel HORSE formalism [18,^|, the radial channel functions 
u-p,-p{ (/cr,Tr) are expanded in terms of the harmonic oscillator functions, 

oo 

%(ro *^^r, '"r) = Yl '^nr(r.) (^r) Rnir (fr) 



(62) 



n=0 



[we omit the indexes J and (±) in notations of the wave functions in oscillator representation: 
anTiVi) (kr) = dntir) i^r) ]■ It is convenient to expand also the internal channel wave func- 
tions ^uirsrixr) in oscillator function series, but generally any other representation for the 
functions ^uirsrixr) can be used. The matrix of the projectile-target interaction within the 
multichannel HORSE formalism is truncated according to eq. (0) in each particular channel 
r independently from the others, and as a result we shall get a set of truncation boundaries 
{A^r}- In the asymptotic region n > N-p, the wave function in oscillator representation is of 
the form 



^(-) 



(+) 



<f(r.)(^r) = 5rr.Q7/(A:r) - Srv^C^J {kr 



(63) 



where cS\k) = Cni{k) ±iSni{k). 

Truncating the Hamiltonian matrix H^^, at n = N^ in each particular channel F and 
diagonalizing the truncated Hamiltonian matrix H^^, , one can calculate matrix elements 



rr' ^ _ \p V7AnJ l\n' p' 
^nn' Z^ /? _ /? n',n'+l i 



A ^A 



(64) 



which we use to construct the matrix 
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/ 



[Q] 






N2 



glM 



Ni Nm 



Q21 Q22 



g2M 



N2NM 



y^NuNi yNMN2 



r.M M 

yNM Nm ) 



(65) 



In eq. ( |6^ ) T^„+i is the relative motion kinetic energy matrix element for the channel F, and 
E\ and 7^^ are eigenvalues and corresponding eigenvectors of the matrix ii^, . Equation ( |6^ 
is a generalization of eq. (|19|) to the multichannel case. The matrix elements Q^i are useful 
for calculations of the multichannel wave function in the harmonic oscillator representation 
in the interaction region n < N^ in any particular channel F by 



anr(r,)(^r) = I]^n^^, «^j„+i,r'(r,) (^rO , (66) 

r' 

while eq. ( |63D should be used for calculation of a^wp^-, (fcr) in the asymptotic regions n > Nj-. 



To make use of eq. (|63|) , one should first calculate the S'-matrix. Within the multichannel 
HORSE formalism it can be done by the following expression [p9[]: 



[s] = {[ci:\k)] - [G] [ci:i{k)]y' {[ci^\k)] - ig] [ci,-i{k)]} 

or, taking into account the symmetry property of the S'-matrixQ, 



(67) 



[S] 



Tr 



(68) 



where the index Tr labels a transposed matrix, 



[S] = [[cVik)] - [d^lAk)] [Gf'} {[Cy\k)\ - [C^lAk)] [GfY' 



(69) 



^The ^-matrix is defined with cuts in the complex energy plane starting from threshold energies 
er of all channels T. If the threshold energies er are different for different channels, the symmetry 
property (^) may not hold in the complex energy plane. In this case some of the formulas 
below should be modified. However, all the formulas below can be used at energies E such that 
Re £' > max {er} or Re i? < min {er}. 
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lal matrices 


c. 


(±) 

V+1 


[k) 


and 


'd^\k) 


are defined as 




c 


fXuAk,) ] 


^N+lW 


= 






C^NXuAk,) 








\ 









^Nm+1,Im (^m) J 



(70) 



and 



/ 



d^\k) 



>(±) 



Ctvi h (^1 



C^nX ik2) 



\ 



.(±) 



(71) 



^v:;/„(m; 

In order to define tlie P-matrix within the multichannel HORSE formalism, we substitute 
the S'-matrix in eq. ( p7b| ) by the right-hand side of eq. (^). After some algebra, we derive: 

[P] = [b] ( \hM {k,b)]' \c'}^\k)] - \h(-) {k,b)]' [cj,+) {k) 



HM{k,b)\ [C^j,~Uk)\-[H(^^{k,b)\ [&^A,{k) )[g] 
x{[HM{k,b)] [Ci^\k)] - [H^-\k,b)] [Ci^\k) 



(72) 



- {[H(+Kk,b)] [ci^Uk)] - [H^-Kk,b)] [dAlik)]) [gf^y . 

P-matrix is known [^ to be a real matrix. Nevertheless, most of the matrices entering ( |72D 
are complex ones. It is useful for applications to rewrite eq. ( [7^ ) in a form that does not 
contains complex variables: 



[P] = [b]^[F{k,b)\ [CN{k)]-[G{k,b)\ [SN{k)] 

- ([Fik,b)]' [CN+iik)] - [Gik,b)]' [SN+iik)]^ [gf 

X {[F{k,b)] [CNm - [Gik,b)] [SNm 

- {[F{k,b)] [C^+i(A:)] - [G{k,b)] [S^^,{k)]) [gf''} 



(73) 



where the matrices 



Fik,b) , G{k,b) , F{k,b) , Gik,b) , [SNik)l [G^m, [SN+im, 



and [GN+i{k)] are diagonal matrices with matrix elements equal to Fi^{kr,br), 
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d 



rr=br 



arr 
and CNr+i,ir{kr), respectively, e.g 



rfrr"'" 



Gi^{kr,rr) 



5'Arpip(fcr), C'Afrtl^r), SNr+l,lri^T)j 



; ^iVptr 



rr=br 



F(fc,6) 



FiAkiM) 



Fi,{k2M) 



\ 







Flj,,i{kM, &Af) / 



(74) 



Eq. ([73| ) generalizes eq. p7| ) to the multichannel case. It can be used for the calculation 
of the multichannel P- or i?-matrix, [R\ = [P]^ , within the harmonic oscillator expansion 
method for any set of the channel radii {6r} large enough to neglect the potential energy at 
the distances r-p > b^ in each channel F. 

P-matrix poles coincide with eigenenergies of the truncated Hamiltonian matrix H^^' 
if the matrix [Q] ^ is not actually present in the inverse matrix in the right-hand side of 
eqs.([7^)-(|73|), i.e. when the condition 



F{k,h) \C^^^{k)\ = G{k,h) [SN+i{k)] 



(75) 



is satisfied. Equation ([75|) is a matrix equation that should be solved to find a set of channel 
radii {^r}- However, all the matrices entering eq. ( ffSD are actually diagonal matrices. Thus, 
eq. (^) reduces to a set of M uncoupled equations ( p9D each corresponding to some channel 
r, and the channel radius 6r for any particular channel F can be found independently from 
the others. As a result, the set of natural channel radii {6°} can be introduced with the 
help of eq. (^) applied to all channels F in turn. Of course, the natural channel radii 6p are 
defined with the same accuracy as in the single-channel case. 

If all of the channel radii 6r are set to be equal to the corresponding natural channel 
radii 6p [or, more generally, if the set of the channel radii {6r} is any solution of eq. ([75|)], 
the general expression for the P-matrix (^) reduces to 



(76) 



[P]= j_[lP\[T]\\G{k,lP)\ [SMm-[nk,h')\ [Cr,{k)] 

' G{k,\P)\\S^^rm - [P(A;,60)l'[Cjv+i(A;)]) [^]| \G{k,\P)WCn^^{k)\ 
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where [T] is a diagonal matrix built up by kinetic energy matrix elements T^ ^ _,_^. 
Eq. ([76|) is a multichannel generalization of eq. (p3|). In the same approximation that was 



—^ sin kvrT ^ and Gi^{kr, rr) 



used to derive eq. (plf), i.e., in the case when one can use approximations Fi^{kY,rr) ~ 

^oosLrr-^) in every channel rat 
the distances r^ > br, eq. (JT^) can be simplified to 

[P] = - [b^] [k] [T] {[Sr^m [Sr^+^m + [C^{k)] [C^+i(A:)] 

[[SN+iik)f + [CN+iik)f) [Gik,W)] [CN+iik)]'' [G] [Cik^b^]]'' [CN+iik)]} - [l] , 

(77) 



where [k] is the diagonal matrix of the channel momenta k-p, and 1 is the unit matrix. 

Further simplification of eq. ( [77D can be performed making use of asymptotic expressions 
( |14bD and ( |15bD of the functions Sni{k) and Cni{k). As a result we get: 



[P] 



Here [ro], 



\n + i)In + i + ^ 



m- 



ro] [Q] [ro. 



+1 



(78) 



{N + l){N + l + l) 



and [(3\ are the diagonal matrices with non-zero ma- 



trix elements equal to rj 



{Np + l){Nr + lT + l) and 



h/{fiYUj) where nr is the reduced mass in the channel F, 



Pt 



'2Nr + It + 7/2' 



1/4 



cos < 2/i;rro 



2 4 



'--In 



(79) 



V2A/'r + /r + 3/2^ 

respectively. A multichannel generalization of eq. (^^ can be derived from eq. (|7^) in the 
limit A^r 



oo: 



[P] 



A^ + - + - 
2 4 



[ro]-' [^] [ro] \b' 



(80) 



A^ + - + - 



is a diagonal matrix with (Nj- + ^ + f ) standing for the non-zero matrix 



where 
elements. 

Equation ( pOf ) defines the multichannel discrete analogue of the P-matrix. Note that if 
the reduced mass /ir is the same in all channels F, the matrix [ro] appears to be proportional 
to the unit matrix. In this case expression (^OD simplifies to 
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[P] 



iVH h - 

2 4 



Qi] _[50] 5[^] J^o]2\ _ jj- 



(81) 



A further simplification of eq. (|8ll) is possible if the channel radius b^ is the same in all 
channels F (for example, if truncation boundaries A^r and orbital angular momenta /r are 
the same in all channels F), because in this case [6°] ^ [Q] [W]'^ = [Q]. 

As in the single-channel case, being derived in the quasi classical limit A^ ^ oo, the 
multichannel discrete analogue of the P-matrix provides a very accurate approximation to 
the exact P-matrix even for very small values of truncation boundaries A^r- This is not 
surprising because only the matrix [Q] entering eqs. (pO| ) and (pT]) mixes the channels, all 
the rest matrixes in eqs. (^0]) and (|81|) are diagonal in the channel space and provide in each 
channel the same accuracy as in the single-channel case. 

The above multichannel P-matrix formalism can be used to derive a multichannel gen- 
eralization of our Coulomb-HORSE method. 

Consider a multichannel scattering in the system with the interaction described by a 
superposition of a short-range nuclear potential V^^^' and Coulomb potential Vpp'"'. Con- 
trary to the Coulomb interaction, the nuclear potential V^Y'^'' couples the channels. The 
wave function is given by eq. (|52D where rji^ = argF(l + /r + ^Cr) and Cr are the Coulomb 
phase shift and Sommerfeld parameter, respectively, in the channel F. In the asymptotic 
region rp > bj-, the radial channel functions u-p,-p{ {kr,rr) are of the form: 



u 



Ji+) 



1 



r(P^)(A;r,rr) = - [G\l> {kr,rr) Srr, - G\;> {kr,rr) ^rr, 



(82) 



where 



G\^\k,r) 



^/I^KC^r) ±^l^^F,iCM) 



(83) 



1 



rwv 



exp 



/ Til ' 

±i Ikr — (\n 2kr \- rji 



and Fi{(, kr) and Gi{(, kr) are the regular and irregular Coulomb function, respectively. 

By analogy with section [IV A| , we introduce an auxiliary short-range potential Vp'^ by 
cutting Coulomb potential Vf°'^'' at distances r^ = &r in the each channel F (it is supposed 
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that 6r > R^^'^^ where R^'^'^^ is the range of the nuclear interaction in the corresponding 
channel) : 



sn 1 ^rr' = V^r'^^^^' + SrrVfr^ ^r < W 






0, 



rr > &r 



6r > Rv'"''^ 



(84) 



In the asymptotic region r-p > br, the radial channel functions u-p,-p\ {kr, rp) obtained with 
the auxiliary potential (^) are given by eq. ( |5^ ) with an auxiliary S'-matrix S^^. standing 
for Srr,- 



u 



r(fcr,rr) = ^ (^i+^fcr,rr) S^ - Hi;\kr,rr) 5rr.) • 



Ji+)Sh 

r(r, 



(85) 



As in the case of single-channel scattering, for the set of channel radii {br} we have an 
equation 



)C 



P 



Sh 



(86) 



where 



jC 



is the P-niatrix of the multichannel Coulomb problem and the P-matrix 



pSh 

corresponds to the auxiliary multichannel problem with short-range interaction (0). The 
coupled channel equations for this auxiliary multichannel problem can be solved within the 
multichannel HORSE formalism, and the corresponding scattering matrix 5''^'', the matrix 



of channel functions 



uJ(+)Sh^j,^r] 



The Coulomb P-matrix 



P 



c 



and the P-matrix P^^, can be calculated. 

can be substituted in eq. (|86D by its expression through 

the multichannel Coulomb S'-matrix [S] given by eq. (|57b| ) with H^"^^ {k, b) replaced by 

diagonal matrices G*^^-* {k^b) defined according to the general rule (|7^). As a result, we 

can derive the following expression for the S'-matrix of the multichannel Coulomb problem: 

-1 

(87) 



[h\ 



[h] 



-1 



pSh 
pSh 



G^+\k,b) 
'G^-\k,b) 



G^+\k,b) 
'G^-\k,b) 



[cf. with the general expression of the multichannel S'-matrix through the P-matrix ( A.l ) 



given in the Appendix]. Equation ( |73D can be used to calculate 



:>Sh 



entering (|57|). 
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It is also possible to derive an expression of the multichannel Coulomb 5'-matrix [5*] in 



terms of the auxiliary S'-matrix 



5' 



Sh 



Using the symmetry of the P-matrix [see eq. ( [A.7| ) 



in the Appendix] we rewrite (p6| ) as 



P 



c 



P 



Sh 



Tr 



1-1 



(88) 



We substitute 



P' 



c 



and 



P 



Sh 



in (^8D by their expressions ( |57b|) through [S] and 



S' 



Sh 



respectively, to obtain 

[b] ([G^~\k,b)\' - [G^+\k,b)]' [S] 
H^~\k,b)] - [S^''] [H^+\k,b) 



G^-\k,b)] - \G^+\k,b)] [S] 



-1 



H^-\k,b) 



S 



Sh 



H^+\kM \\b] 



r^-l 



(89) 



Solving equation (|^) with respect to [S] we derive: 



[S] 



r/i^-l 



X 



(90) 



[h] 



-1 



[b] N [H^-\k,b)\ [G^-\k,b)\ - [H^-\k,b)\ [G^-\k,b) 

H'^+\k,b)] [G'^-\k,b)]' - [i^(+)(A;,6)]'[G(-)(A;,6)]| [5^^" 

i^(")(A;,6)] [G^+\k,b)]' - [H^-'^{k,b)]' [G^+\k,b) 

H^+\k,b)] [gW(A;,6)]'- [H'^+\k,b)]' [G^+\k,b)]\[S^'' 
In the derivation we made use of the symmetry of the S-matrix ( pBD (note the footnote |I| 
on page P7|). 

Equation (^) is a multichannel generalization of eq. (|48|) . 

As in the single-channel case [cf. eq. (^)], to obtain the radial channel functions 
MpL{ (kr,rr) of the multichannel Coulomb problem in the interaction region rr < br, one 
should renormalize the radial channel functions Wp^p ) (A;p,rr) obtained with the auxiliary 

uJi+) (^/j^ f^ in the interaction region 



potential (0). The matrix of radial channel functions 
can be calculated as 



u^(+) (k, r)] = [Af] 



u 



J{+)Sh 



{k,r) 



(91) 



where 



^j{+)Sh (j^^ ^-j jg ^YiQ radial channel function matrix of the auxiliary multichannel 



problem, and the renormalization matrix [A/] is given by 
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[AT] = { [gW {k, b)] [S] - [g^-^ {k, b)] } { [i^(+) {k, b) 



S 



Sh 



H^-^{k,b)\} 



(92) 



Note that the renormahzation matrix [Af] is not diagonal. This is because the long-range 
Coulomb interaction modifies in a different manner radial channel functions in different 
channels. As a result, the relative weights of channel amplitudes in the interaction region 
are changed by the Coulomb interaction. This requires some rearrangement of the auxiliary 
wave function in the interaction region that can be achieved only by mixing its components 



u 



J{+)Sh 

r(r,) 



(A;r,Tr), which is generated by non-diagonal matrix elements of [Af]. 



Ji+) 



In the asymptotic region rp > &r, the radial channel functions Wpjp ) (A;r,''"r) are given 



by expression (|82D . 

We suppose to present the results of applications of the multichannel Coulomb-HORSE 
formalism to nuclear reactions in future publications. 



VI. SUMMARY 



In this paper we have developed P- and i?-matrix formalism within harmonic oscil- 
lator representation of scattering theory. With the help of this formalism, one can use 
well-developed P- or /^-matrix technique HH] in combination with standard variational ap- 



proaches based on harmonic oscillator expansion |^0[. In particular, one can use the results 
of nuclear structure calculations obtained by standard shell model codes to find resonance 
positions and widths, cross sections of various nuclear reactions and other scattering observ- 
ables. 

We derived expressions for calculation of R- or P-matrix at any distance from the origin 
covered by the set of oscillator basis functions employed in the variational procedure. We 
have shown that there is a natural channel radius associated with the given oscillator basis. 
For the natural channel radius, the expressions for the P- or i?-matrix are simplified and a 
discrete analogue of the P-matrix can be introduced. The discrete analogue of the P-matrix 
is expressed through the finite difference of the wave function in the oscillator representation 
in the same manner as the usual P-matrix is expressed through the derivative of the wave 
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function in the coordinate representation. It has the same properties as the conventional 
P-matrix, in particular, its poles are just the eigenenergies obtained in variational approach 
with oscillator basis. The discrete analogue was shown to be equal to the P-matrix in 
the quasiclassical limit of large number of basis functions, however, it gives a very good 
approximation for the P-matrix even in the case when only one oscillator basis function is 
involved in the variational procedure. 

The discrete analogue of the P-matrix is useful for calculation of scattering and reac- 
tions of uncharged particles. It cannot be used for scattering calculations in the case when 
Coulomb or another long-range interaction is present in the open channel. However, in this 
case one can calculate the P- or P-matrix at shorter distances than the natural channel 
radius for investigation of scattering observables within the variational oscillator-basis ap- 
proach. For such investigations, we developed a very efficient Coulomb-HORSE formalism. 

In Section 0, we generalized all our results on the case of arbitrary number of channel. 
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APPENDIX: SYMMETRY OF MULTICHANNEL P-MATRIX. 



The expression of the multichannel P-matrix through the S'-matrix is given by eq. ( p7b|) 
which can be used to express the S-matrix through the P-matrix: 



[S] = Ubr^lP] [H(+\k,b)] - [H(+\k,b)]'\ 
[b]-^ [P] [H(~^{k,b)] - [H(-^k,b)]'\ 



(A.l) 
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Using the symmetry of the 5'-matrix (^) (note the footnote |T| on page ^^ , we can rewrite 
eq. ([A. 1|) as 

' "i7(-)(A;,6)] [Pf' [h]-^ - [H^-\k,b) 
'H^+\k,b)] [Pf'' [b]-' - [H^+\k,b) 
It follows from eqs. ( |A.1| )-( |X^ ) that 



[S] 



(A.2) 



[b]-' [P] \H^+\k,b)] - \H^+\k,b)]'\ \[b]-' [P] \H^-\k,b)] - \H^~\k,b) 



H^-\k,b) [P] W - H^-\k,b) H^+\k,b) [Py^'ibr^- H^+\k,b) 



(A.3) 



(A.4) 



and after some algebra we obtain: 

'H^+\k,b)]' [H^-\k,b)] - [H'^-\k,b)]'[H'^+\k,b)]\[Pf''[b]-^ 

= [b]-^ [P]l[H^+\k,b)]' [H^-\k,b)] - [H^-\k,b)]' [H^+\k,b) 
The Wronskian of the solutions H} (k, r) and Hi {k, r) defined by eq. (^5|), is 

Hi+^\k,r) Hl~\k,r) - Hj-^' {k,r) Hl^\k,r) = ^ ^^ . (A.5) 

Therefore eq. (|A.4|) appears to be equivalent to 



[pr = [p] 



(A.6) 



or 



[P] 



[P] 



Tr 



(A.7) 



where the diagonal matrix 

/xr- 

It is clear that the multichannel P-matrix is symmetric 



[ji] [6] and [ji] is the diagonal matrix of reduced masses 



[p] = [py 



(A.8) 



if 



is proportional to the unit matrix, i.e. if the ratio -— is the same in all channels F. 

br 
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FIGURES 
FIG. 1. Relative deviation (60 — b)/b of the natural channel radius 60 from the exact solution 

b of eq. (^) vs E for various values of N . Solid lines — s waves; dash lines — p waves; dot-dash 

lines — g waves, fioj = 18 MeV, reduced mass m corresponds to the neutron scattering by ^ = 15 

nucleus. 

FIG. 2. P-matrix at the natural channel radius for neutron scattered by A = 15 nucleus vs 
energy in the cases A^ = 1 (a) and A^ = 9 (b). Solid line — exact P-matrix calculated by (P7|); dash 
and dot lines — discrete analogues of the P-matrix calculated by ( p8| ) and (^5[), respectively; up 
and down arrows indicate positions of the poles of the exact P-matrix and its discrete analogue, 
respectively. Calculations are performed with huj = 18 MeV and Woods-Saxon potential of ref. [^. 

FIG. 3. s-wave phase shift 5q of the p-^^N scattering as a function of the channel radius h used 
for construction of the potential V^^ at energies E = 2 MeV (a) and i? = 10 MeV (b). The hori- 
zontal solid line correspond to the exact value of 5q. Phase shifts calculated with fiuj = 18 MeV and 
A^ = 4, A^ = 9 and A^ = 19 are presented by short-dash, solid and long-dash curves, respectively; 
dot-dash (dots) curve presents the results obtained with fiuj = 10 MeV and A^ = 9 {huo = 26 MeV 
and A^ = 9). 

FIG. 4. (a) s-wave phase shift Sq of the p~^^N scattering as a function of energy E. Solid curve 
shows the exact values of 6q. Calculations by the method suggested in this paper with hu) = 18 MeV 
and channel radius 6 = 7 fm, and A^ = 10, A^ = 8, A^ = 6 and A^ = 4 are presented by the long-dash, 
dot-dash, dot-dot-dash and short-dash curves, respectively. Calculations by the method suggested 
by Kiev group [^ are presented by the thin solid curve with diamonds (A^ = 50 and A^ = 70) 
and by dots (A^ = 8 and A^ = 70). (b) The same but in the vicinity of the resonance only. 

FIG. 5. (a) Wave function in the oscillator representation squared, a^^i^k), as a function of 
energy for Z=0 and n = 0, 1, 2. (b) The same but in the vicinity of the resonance only. See fig. ^ 
for details. 
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FIG. 6. Radial wave function ui{kr) for s-wave p-^^N scattering at energies E = 3 MeV (a) 
and E = 15 MeV (b). Solid curve — exact; calculations by the method suggested in this paper 
with Tllo = 18 MeV and channel radius 6 = 7 fni: dash curve — reconstruction with N = M = 10; 
dots — reconstruction with A^ = 10 and M = 100 [see eq. ([5l|)]. 

FIG. 7. p-wave phase shifts of the p-^^N scattering (a) and the corresponding wave functions 
in the oscillator representation squared, a^;(A;), for 1=1 and n = 1 (b), as functions of energy. See 
fig. § for details. 

FIG. 8. The same as fig. but for d-waves. 
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